EXPONENTIAL MOMENTS OF AFFINE PROCESSES 



MARTIN KELLER-RESSEL AND EBERHARD MAYERHOFER 

Abstract. We investigate the maximal domain of the moment generating 
function of afiinc processes in the sense of Dufiie, Filipovic and Schachermayer 
■ [2003] and we show the validity of the affine transform formula that connects 

exponential moments with the solution of a generalize d Riccati differential 
equat i on. Our result extends and unifies preceding ones bvlGIasse rman and KimI 
I2OI0II. IF ilipovic and M averhofe'rl l2009l l. lKallsen and Muhle-Karbd I2OI0II and 
[Spreii and VeermanI |2010t in that it allows processes with completely general 
jump behavior, applies to any convex state space and provides both sufficient 
and necessary conditions for finiteness of exponential moments. 
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1. Introduction 



This article investigates the maximal domain of the moment generating function 
of an afBne process. An affine process is a time- homogeneous Markov processes 
X on a finite-dimensional state space D C M'' whose characteristic function has 
the following property: There exist a complex-valued function (j> and a C'-valued 
function ip such that 

(1.1) <i>{t,u,x) :=E[e<"'-^'> | Xq = x] = e-^^'^^^+W'^'-^^'^^, 

for all u G iM.'^, t > and x <E D. This so-called affine property implies that the 
Forward-Kolmogorov-PDE for $(t, u, x) 

d 

— u, x) = A^{t, u, x), $(0, u, x) — exp((w, a;}), 

where A denotes the infinitesimal generator of X, can be reduced to a system of 
non-linear ODEs, commonly referred to as generalized Riccati differential equations, 
which are of the form 



(1.2a) 
(1.2b) 



^^cj){t,u)^F{^P{t,u)), </)(0,w) 



0, 



R{'ijj(t,u)), ■)jj{0,u) = u. 



A natural and important question is whether formula (jl.ip and the generalized 
Riccati system (|1.2p can be extended to real exponential moments {u G M.'^) or 
complex exponential moments {u G C*). One might expect that if F and R can 
be suitably extended - e.g. by analytic extension - then the exponential moment 
|^g('":-^T>j jg finite if and only if a solution to the extended Riccati system exists 
up to time T, and that in this case also (jl.ip remains valid. A statement of this type 
is usually referred to as affine transform formula. Showing such a formula in full 
generality is far from trivial - difficulties include the fact that analytic extension 
of F and R may not be possible, that solutions of the extended Riccati equations 
might not be unique, and that the differentiability of i h-s- 4>{t,u) and t i-^ tp{t,u) 
is not obvious from (|l.ip . The latter problem of showing that differentiability of 
(/) and tp can be concluded from the definiti on of an affine processes is known as 
the re g ularity problem for affine processes (cf. Duffie. Filipovic. and Schachermaver 



2003l |. lKeller-Ressel. Teichmann. and Schachermaveil f2011b| . lcuchierol 120111 ) 

Several articles have been concerned with showing the affine transform formula 
under different conditions on the process X or the state space D. In particular we 
mention the following contributions: 

Glasserman and Kim 2010| | show the affine transform formula for affine 
diffusion processes on D = R~^n x I^" under a mean-reversion condition; 
Filipovic and Maverhofer |2069t show the affine transform formula for affine 
diffusion processes on D = x K"; 

Kallsen and Muhle-Karbe 2010l | show that for affine semi-martingales on 



a general state space existence of a solution to the extended Riccati sys- 
tem implies finiteness of exponential moments under very mild additional 
conditions; 

Spreij and VeermanI 2010l | show the affine transform formula for affine pro- 



cesses whose jump measures possess exponential moments of all orders, and 
for arbitrary convex state space; 
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Cheridito and Wugalteil j2012l | show the afhne transform formula with kiUing 
for affine processes whose jump measures possess exponential moments of 
all orders. 

In this article we generalize and unify most of these results by showing the affine 
transform formula for real exponential moments of an affine processes on an arbi- 
trary convex state space, without imposing any conditions on the jump measure. 
In addition we show the affine transform formula for complex exponential moments 
for affine processes on D = x R" (the canonical state spaces) and on D = St 



the p ositi ve semidefinite d x d matrices). These two state spaces (see iDuffie et al 
2003| and Cuchiero. Filipovic. Maverhofer. and Teichmann 201lj ) are of particu- 



lar interest both from the theoretic viewpoint and from the applied one. 

The outline of this paper is as follows: In Scction[2]we present general definitions, 
some useful notation and our main results: 

• Theorem 12 . 141 proves the affine transform formula in terms of minimal solu- 
tions to the so-called extended Riccati system, which comes from considering 
()1.2p in the real domain. Here we only require the state space to be closed, 
convex and with non-empty interior. The proof of Theorem l2.14l is provided 
in Section m 

• Theorem 12.261 extends the validity of the affine property (jl.ip to complex 
moments u — p + iz, where z € M.'^. This extension succeeds under the 
premise that the p-th real moment is finite, or equivalently, that the ex- 
tended Riccati equations are solvable until time T. The result holds for the 
state space M™q x R" and - under some mild additional conditions - for 
the state space S'^. For the proof of Theorem 12.261 see Section [5] 

In Section [3] several applications of our results to mathematical finance are out- 
lined. Finally, Sections |4] and [5] contain the proofs of our main results for real 
moments and complex moments respectively. 

2. Definitions and Main Results 

2.1. AfRne Processes. Let (17, F) be a filtered space, with F = {J-t)t>o a right- 
continuous filtration. We endow R'', (d > 1) with an inner product (.,.) and let 
D be a non-empty convex subset of R'', which will act as the state space of the 
stochastic process X we are about to define. The state space D has a measurable 
structure given by its Borel a-algebra B{D), and without loss of generality - see the 
explanation after Definition 12.21 - we may assume that D contains and that the 
linear span of D is the full space R'^. Under this assumption it follows in particular 
that the interior D° of D is non-empty. Associated to D is the set 

(2.1) U = ^ueC^ -.x^ e^'''"^ is bounded on d| . 

Finally let {P'^)xeD be a family of probability measures on the filtered space 
{il,J-',¥) and assume that T is P^'-complete for every x £ D. 

Let A be a cadlag F-adapted time-homogeneous conservative Markov process with 
state space D. More precisely, writing 

(2.2) pt(a;,A) =P^(At e A), (t > 0,x e D, A e B{D)) 
for the transition kernel of A, pt{x, A) satisfies the following: 

(a) X i-> pt(x. A) is ;B(£>)-measurable for all t > 0,A £ B{A), 
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(b) pt{x, D) = 1 for all t > 0, .t G Z), 

(c) Pi){x, {x}) = 1 for all x G D, and 

(d) the Chapman-Kolmogorov equation 



Pt+s ix,A) = / ptiy, A)ps {x, dy) 
Jd 

holds for every t,s >0 and (x, A) e D x B{D). 



Remark 2.1. Since X is cadlag, the law of X under is a probability measure on 
the Skorokhod space of cadlag paths ©(R^o, ^'^), for each x G D. There will be no 
loss of generality by directly interpreting as a measure on this path space. 

Definition 2.2 (AfFine Process). The process X is called affine with state space 
£), if its transition kernel pt{x,A) satisfies the following: 

(i) it is stochastically continuous, i.e. lims_j.t .) = pt{x,.) weakly for all 
t>0,xGD, and 

(ii) there exist functions cf) : M^o x — > C and -0 : M^o xU such that 
(2.3) / e<"'^>pt(x, dO = exp {<j>{t, u) + {x, u))) 



JD 

for all t > 0,x E D and u eU. 

Remark 2.3. We explain why it is no loss of generality to assume that D contains 
and linearly spans the whole space W^: For an arbitrary non-empty convex subset 
D of K."^, let afF(D) be the smallest affine subspace of that contains Z?, and let 
(xo, xi, . . . , Xk) be an affine basis of aff(D) such that xq g D. Let h : a.S{d) R*^ : 
X I— > A^ {x — xq) be the projection to canonical affine coordinates, i.e. h{xo) = 
and h{xi) = e, for each j e {1, . . . k}. Set D ^ h{D) C K*^ and X = h{X). Then D 
is convex, contains and linearly spans R*^ . It is easily verified that X is again an 
affine process with 

(2.4) 4)[t, u) = (j){t, Au) + {xo,ip{t, Au) - u) 

(2.5) i}{t,u) = A+i}{t,Au), 

where A'^ is the Pscudoinvcrse of A (or any other k x c?-matrix such that A^ A = 
idfc). 

The next result shows that an affine process is a semimartingale with affine 
(differential) semimartingale characteristics. 



Theorem 2.4 (|Cuchierd 2011|). Let X be an affine process with state space 
D C R''. Then, for each x G D the process X is a -semimartingale with semi- 
martingale characteristics 

(2.6a) At= f a{X,_)ds, 



(2.6b) 



Bt= f b{Xs-)ds, 

JQ 



(2.6c) i'{uj,dt,d£,) ^ K{Xt-{uj),d^)dt 
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where a{x),b{x) and K{x, d^) are ajfine functions of the form 

(2.7a) a{x) — a + xia^ + • ■ • + Xdct'^, 

(2.7b) b{x)^b + xil3^ + --- + XdP'^, 

(2.7c) K{x, d^ ^ m{d^) + xin\dO + ■■■ + a;d/(dC) 

and for each x ^ D it holds that a{x) is a positive semidefinite d x d matrix, b{x) 
is a -vector, and K{x,d^) is a Radon measure on , satisfying 



f (lief Al) A-(x,dC) <oo 



and K{x,{0}) = 0. 



Proof Follows from ICuchierol |201lL Thm. 1.4.8 and 1.5.4]. □ 



Remark 2.5. Note that several of the assump tions mad e at th e beginning of the 
section could be slightly weakened: Following Cuchiero |2011 any affinc process 



has a cadlag modification; moreover the P^-completion of the filtration generated 
by an affine process is automatically right-continuous. 

2.2. Real Moments of AfRne Processes. 

Definition 2.6. Given an affine process X and the associated functions 
{a{x),b{x),K{x,d£^)) in (|2.6p . define for each x £ D the function TZ^ : M'' — )■ 

(— oo, oo] by 

(2.8) TZ,iy)^Uy,aix)y) + {bix),y)+ [ (e<«'^^> - 1 - (/i(e), y)) Kix^d^, 



where hi^) = '^{\(\<i}i- 



For each fixed x G D, the function TZx is a convex and lower semicontinuous 
functior{3 that may take the value +oo. As for any convex function, the effective 
domain is the set of arguments for which TZ^ takes finite values. Taking the 
intersection over all a; G I? leads to the following definition. 

Definition 2.7. Given an affine process X and the associated function TZx as in 
Def. define 



(2.9) y= Pi {yem.'^: / e'^y^^'>K{x, dO < oo 



As an intersection of convex sets, also y is convex. Moreover, y contains and 
hence is non-empty, because TZx{0) — for all x £ D. 

Since the functions a(x), b{x) and K{x, d^) are affine in x, we can decompose TZx 
into TZx{y) = F{y) + {R{y),x). For arguments y & y, the functions F and R are 
uniquely specified, since D contains and d linearly independent points. 

Proposition 2.8. Let X be an affine process with state space D. Then there exist 
functions F : y ^ R, R : y ^ such that 

TZx{y)=F{y) + {R{y),x) 



"'^lower semicontinuity follows from Fatou's Lemma applied to the integral with respect to K{x, d^). 
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for all X £ D, y E y. Let (ei, . . . , erf) be a set of d linearly independent vectors in 
. Then we can write F and Ri{y) := {R{y),ei) as 

(2.10a) Fiy) = i {u, ay) + {b, y) + [ (e<«-^> ~ 1 - {hiO,y)) m(dO , 

(2.10b) R,{y) = i (y, a'y) + {p\y) + / (e<«-^> - 1 - {h{i),y)) y^\d£,), 
^ Je<'\{o} ^ ' 

with /i(0 = l{|5|<i}C 

Proof. Follows immediately from Definition 12.61 and Theorem l2.4l □ 

Remark 2.9. Setting x = in (|2.8p yields that F{y) is a convex and lower semi- 
continuous function of Levy-Khintchine form. The same is not necessarily true for 
. . . , Rd, since the matrices a* may not be positive semidefinite or the measures 
/i* may be signed measures. 

We use the functions F{y) and R{y) to set up a system of ODEs associated to 
the affine process X. These equations play a key role in our main result. 

Definition 2.10 (Extended Riccati system). Let X be an affine process and F,R 
and y be defined as in Definition 12. 71 and Proposition l2.8l Let T > Q,y € y and let 

P-t^p{t,y), q:th^q{t,y) 

be C^-functions mapping [0,r] to K resp. 3^ that satisfy 

(2.11a) ^^p{t,y)^F{q{t,y)), p{Q,y) = Q, 

(2.11b) ^^q{t,y)^R{q{t,y)), qiO,y)^y 

for all t G [0,T]. Then we call {p,q) a solution (up to time T and with starting 
point y) of the extended Riccati system associated to X. 

It is important to note that in general the function R is locally Lipschitz con- 
tinuous only on the interior of y, but may fail to be Lipschitz continuous at the 
boundary of y. Hence solutions of (|2.1ip reaching or starting at the boundary of 
y may not be unique. For this reason we add the following definition. 

Definition 2.11 (Minimal Solution). Let X be an affine process, and let {p,q) 
a solution up of T starting at ?; G 3^ to the associated extended Riccati system. 
We call (p, q) a minimal solution, if for any other solution (p, g) up to T < T and 
starting at the same point q{0, y) = g(0, y) = y it holds that 

(2.12) p{t, y) + {q{t, y), x) < p{t, y) + {q{t, y),x) 
for aU t e [0, f] and x e D. 

Remark 2.12. By setting (72;(t, y) := p{t,y) + (q(t,y),x)., the extended Riccati system 
may be written in condensed form as 

(2.13) ^^qAt,y)^n^{q{t,y)), g,(0,y)=y, Wx G D. 
In this notation the minimality property can we written as 

qAt,y)<qAt,y), VxeD,te[o,f],yey. 
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Remark 2.13. The following properties are easy to see: If for a given starting 
value y G y there is only one solution to the extended Riccati system, then it is 
automatically a minimal solution. Also, if for a given starting value a minimal 
solution (p, q) exists up to time T, it is automatically the unique minimal solution. 
Indeed, if there were another minimal solution (p, q) , then 

P{t, y) + y), x) = p{t, y) + {q{t, y),x) 

for all t G [0, T], x G D. Since D contains d linearly independent points and it 
follows that p = p and q ~ q in this case. 

We can now formulate our main results on the behavior of exponential moments 
of affinc processes. 

Theorem 2.14 (Real moments of Afhne Processes). Let X be an affine process on 
D and let T > 0. 

(a) Let y G Mf^ and suppose that [e<y^^^>] < oo for some x e D° . Then y ey 
and there exist a unique minimal solution {p, q) up to time T of the extended 
Riccati system (|2.1ip . such that 

(2.14) le'^y'^'A = exp {p{t, y) + {q{t, y), x)) 

holds for allxe D,t(E [0,T]. 

(b) Let y € y and suppose that the extended Riccati system (|2.1ip has solutions 
{p, q) that start at y and exist up to T. Then E^ [e^*''"''"^^] < oo and there exist 
unique minimal solutions (p, q) up to time T of the extended Riccati system 
such that (Pl^ holds for allx e D,te [0,T]. 

Remark 2.15. We emphasize that in point (b) of the theorem p = p and q = g does 
not necessarily hold, i.e. the candidate solutions (p, q) have to be replaced by the 
minimal solutions {p, q) in order for (|2.14p to hold true. 

The following Corollary is a conditional version of Theorem 12.141 

Corollary 2.16. Suppose that the conditions of either Theorem \2.1J^ or \2.l]j^ are 

satisfied, and let {p, q) be the associated minimal solutions of the Riccati system 
((2ll|) . Then also [e(9(T-t,y),^t>] < oo and 



E^ 



c^^{p{T -t,y) + {q[T ~t,y,Xt)) 



holds for all x e D,t(E [0,T]. 



The next proposition provides a way to identify whether some solution (p, q) of 
the extended Riccati system is in fact the minimal solution. 

Proposition 2.17. Let X be an affine process and let {p,q) be a solution up to 
time T > of the extended Riccati system associated to X . Each of the following 
conditions is sufficient for {p, q) to be the unique minimal solution: 

(a) X is a diffusion process; 

(b) y = W; 

(c) y is open; 

(d) qit,y)ey° for allte[0,T]. 
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Proof. From Definition 12.71 of y it foliows tliat (a) (b) (c) => (d), i.e. it is 
sufficient to sliow tliat (d) implies uniqueness of tlie solution (p, q). But R is locally 
Lipschitz on y° , such that standard ODE results imply that {p, q) is the unique 
(and hence unique minimal) solution of the extended Riccati system (|2.1ip . □ 

Remark 2.18. Condition (b) is equivalent to /|j|>i e^^'^'^ K{x, d^) < oo for all x £ D, 
y G R'', which is often referred to as the jump measure having exponential moments 
of all orders. In thi s spec ial c a.se analogues of Theorem 12.141 have been shown in 
Spreii and Veerman 20ld | and Cheridito and Wugalter 2012| . 



We briefly discuss two important special cases, in which great simplifications of 
the results occur. These cases have been treated previously in the literature, but 
serve as a first 'sanity-check' of the main results of this article. 

Example 2.19 (Affine Diffusion). Suppose that the affine process X is a diffusion. 
In this case K{x, .) = ior all x £ D and consequently 3^ = R'' and the functions 
F{y), Ri{y), . . . , Rdiy) are quadratic polynomials (hence locally Lipschitz continu- 
ous everywhere). In this case any solution of the extended Riccati system is unique 
and there is no need to introduce the concept of minimal solutions - see Propo- 
sition 12.1161 above. Thus, Theorem 12.141 holds true even with 'minimal solution' 
replaced by 'solution'. For the case of affine diff usions on canonical state spaces , 
the analogue of Theorem 12.141 has been shown in Filipovic and Maverhofer 20091 
Theorem 3.3]. 

Example 2.20 (Levy process). Suppose that X is a Levy process. Then X is 
an affine process with R{y) — and with F{y) equal to the Levy exponent of 
X. Consequently y is simply the effective domain of the Levy exponent. The 
extended Riccati system has unique global solutions for each y G y, which are 
given by p{t,y) = tF{y) and q{t,y) = y ioi t > 0. It follows from Theorem 12.141 
that [e^^'^^')] is finite if and only if y G 3^, and in case of finiteness we have 
|^e(a.^'t>j — cxp{tF{y) + {y,x)). In particular, finiteness of exponential moments 
is a time-independent property, i.e. for given y G R'' the exponential moment 
grr ^e'^v,Xt)-^ jg gjjji-j, fQj. aU t > or for no i > 0. Of course, all these results 

are w ell-known in the case of Levy processes and can be found for example in lSatol 
19991 Thm. 25.17]. 



2.3. Complex Moments of Affine Processes. In this subsection we give an 
analogue of Theorem 12 . 141 for complex exponential moments of X. The first step is 
to analytically extend the functions F and R. We introduce the following notation: 
For a set A C R'' write 

S{A) :^ {u e : Rc u e A} 
for the complex 'strip' generated by A. 

Proposition 2.21. Let X be an affine process and suppose that y° ^ 0. Then, for 
every x E D, the function TZx defined in (|2.8p has an analytic extension to S{y°) 
which we also denote by TZ^ ■ Moreover it holds that 

TZx{u) = F{u) + {R{u),x) , xeD,ue S{y°) 

where F, R are the analytic extensions of the functions defined in p.lOj) to S{y°). 

Proof . Follo ws from standard results on Levy-Khintchine type functions; see, e.g., 
SatoL Il999l Theorem 25.17] □ 
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Definition 2.22 (Complex Riccati System). Let X be an affine process such that 
y° ^% and let F, R be defined as in Prop. [2311 Let T > 0, u G 5(3^°) and let 

(j) : t 1^ (j}{t, u), ip : t 1-^ ip{t, u) 

be C^-fmictions mapping [0,r] to C resp. 5(3^°) that satisfy 

(2.15a) ^^<j>{t,y) = Fm,y)), 0(0, y) = 0, 

(2.15b) ^7^(t, y) = Riiit, y)), V(0, y) ^ y 

for all t G [0, T]. Then we call {(t),ip) a solution (up to time T and with starting 
point u) of the complex Riccati system associated to X . 

Remark 2.23. Let us compare the complex Riccati system to the extended Ric- 
cati system (|2.11ap - ()2.11bp . We observe that if w G S{y°) is real valued, i.e. has 
Re u = y and Im u = 0, then any solution (0, tp) up to time T of the complex Riccati 
system is also a solution of the extended Riccati system, i.e. setting p{t, y) — (j){t, u) 
and q{t, y) — ■0(t, u) for all t G [0, T] defines a solution (p, q) of the extended Riccati 
system. The reverse is not necessarily true. Furthermore we point out that for a 
given starting value u any solution (0, -0) of the complex Riccati system is auto- 
matically the unique solution. This is in contrast to the extended Riccati system, 
where solutions starting at the boundary may be non-unique. This difference is just 
a consequence of the fact that solutions of the complex Riccati system are restricted 
to stay in the open domain 5(3^°), on which F and R are locally Lipschitz. 

Assumption 2.24. Let X be an affine process with state space D and assume 
that either 

(i) D = X E", or 

(ii) D = and there exists some x G S^^ such that a{x) either vanishes, or it 
is non-degenerate. 

Remark 2.25. Note that in the notation of p.7|) a{x) is given as a symmetric 

d{d+lX ^ d{d+l]_ Qf 

course we can also interpret it as quadratic form on 
which is more natural and, in particular, a coordinate free notion. A simple 
characterization of du]) in terms of the admissible parameter set is given in Remark 

Em 

The analogue of Theorem 12.141 for complex moments reads as follows. 

Theorem 2.26 (Complex Moments of AfSne Processes). Let X be an affine process 
that satisfies Assumption 2.24\ Let T > 0, u £ S{y°) and suppose that the extended 



Riccati system (j2.11ap - (j2.11b|) has a solution (p, q) with initial value Re u up to time 
T such that q{t,Kcu) G 3^° for all t G [0,r]. Then also the complex Riccati system 
p.l5p has a solution (0, ip) with initial value u up to time T, [|e<"^-^'>|] 
and 

(2.16) [e<"-^'>l = exp (0(i, u) + {i){t, w), x)) 



for allx e D,te [0,T]. 
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3. Applications in Mathematical Finance 



This section presents applications of our main r esults, Theoreml2.14l and Theo- 
rem [51111 to mathematical finance in the spirit of Duffie et all l2003l . Section 13]. 
We consider the following generic setup: A traded asset S is modeled by the expo- 
nential of an affine factor process X with state space D, i.e. S = e<^'^> for some 
6 G M."^. Moreover, bond prices are given through an affine short rate model of the 
form 

rt^L{Xt) = l + {X,Xt), 
where I £ M and A £ R*^. This setup includes in particu l ar affine term st r uctur e 



models of interest rates ( Cox. Ingersoll. and Rosd jl98 



Dai and Siiig^ 
Bated |2000l |r 



eton 



.^.^^ ...^v...^... ^..^ ..-.^^ j^^^j, .Duffie and Kan 1996 1. 

2OO0I, etc.), affine stochas t ic vol atilitv models (jHeston 1993 1, 



Barndorff-Nielsen and Shephardl |2001] . etc.) and combinations with 



possible correlation of short rate and asset prices. Also credit risk can be included, 
when rt is interpreted as a superposition of a risk-free short rate and an affine 
default intensity process (cf. iLando 1998l |). Moreover, we can cover a setup with 
multiple possibly dependent assets simply by setting S" = exp {Oi,X) for different 
9i G Mf^. For most applications the measures {V^)x£d should be considered risk- 
neutral measures, although there are few cases where also the behaviour under the 
physical measure is of relevance. Many problems of interest can be reduced to 
determining the J^t-conditional expectations 



(3.1) 



iT.t 



g(x)=E-[e-/"^(^=)'^^5(AT) I -Ft], 



for some measurable function g : D ^ R. In particular 

• (7=1 corresponds to bond pricing; 

• g{x) = e^^'^^ corresponds to checking for the martingale property of the 
discounted asset price; 



y g corresponds to calculating expectations of the type 



[S^ \ which are relevant for evaluation of power utility and determining 
the time of 'moment explosions'. 



^{u,x) 



u £ C corresponds to Fourier methods for the pricing of 



European contingent claims. 

For a more detaile d account of the lit erature on affine processes in financial mathe- 
matics, we refer to lDuffie et al.l 2003L Section 13]; for an easy-to-read introduction 



to dis counting and pricing techniques (using the Fourier-Laplace transform), we re- 
fer to iFilipqyic and Maver hofed j2009l Section 4]. Let us also remark that already 
Duffie et all |2003l Section 11] gives sufiicient conditions on an affine process such 
that the pricing operator Qt,* is well defined, but the results only apply to the 
state space M™q x R" and conditions are less general than the ones we obtain. 

To deal with the di scounting te r m in (|3.ip we use the extension-of-statc-space 
approach outlined in [Duffie et al.l . l2003l Section 11.2]. We define the extended 
state spac^ D := D x R. Let (a, a, b, P, m(dd. ^(d^)) be the parameters of X in 



the sense of Theorem 12.41 Following jPuffie et al.l . l2003l Section 11.2] we have that 



Z := {X, Y) where Yt :~ y + /q L{Xs)ds is an affine process on D with parameters 



alternatively, one could use the minimal spaee D X Img(L), where Img denotes the image of D 
under the map L. 
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{a',a',b',|3',m'{dO,^^'{dO) given by 



a 






( a. 


M 






\ 







K 


j 







and 

* = l,...,d, /3i+i = 0, 

and finally 

where ^0(0^?') denotes the unit mass at 0. Let F{u) and R{u) be the functions 
associated with X through Proposition 12.81 and let q G C. Then we can introduce 
the new functions 

F'{u, q) F{u) + Iq, R'{u, q) = R{u) + \q, 

which are related to the functions (F^, Rz) of the extended process Z in the way 
that Rz = {R'{u,q),0), while Fz{u,q) = F'{u,q). We consider now solutions 
(p{t, u, q) and ijj{t, u, q) of the system 

(3.2a) dt^{t,u,q) = F'{l^{t,u,q),q), (j){0,u,q) = 0, 

(3.2b) dtip{t,u,q) = R'{tp{t,u,q),q), -0(0, w, g) = u. 

Note that ip still is d-dimensional. These solutions are related to the (not necessarily 
unique) solutions 0^, ipz of the corresponding {d+2) dimensional system associated 
with Fz, Rz as follows: (jjzit, (u, q)) = (t>{t, u, q) and ipzit, (u, q)) = {ipit, u, q), q). 

3.1. Bond Pricing in AfRne Term Structure Models. The follo wing result is 



an im mediate consequence of Theorem 12.141 As such it gener alizes Duffie et al. 



2003 . Proposition 11.2] as well [Filipovic and Maverhofer . I2OO9I Theorem 4.11 



Theorem 3.1. Let r > 0. The following are equivalent: 

(1) W[e- ^0 ^'^'>)<^'] < 00, for some x e D° . 

(2) For q = —1, there exists a solution {4>,')p) on [0, r] to the generalized Riccati 
differential equations p.2ap - (j3.2bp with initial data it = 0. 

In any of the above cases, let us define A(t) :— ~4>{t, (0, —1)), B{t) := —tp{t, (0, —1)) 
from the unique minimal solution {4>,ip) of equations p.2ap - ()3.2bp FI Then the price 
P{t, T) of a zero-coupon bond is given, for all < t < T < t, and all x £ D, by 

(3.3) Pit,T) := E-[e- J''^^^)''^ | Ft] = e-^(^-*)-<s(^-*) A) 



3.2. Martingale Conditions. Conditions for the exponentials of affine processes 
to be martingales have been obtained e.g. in lMaverhofer. Muhle-Karbe. and Smirnov 
(2nii |. The following result extends known criteria and follows again from Theo- 
rem UTTH 

Theorem 3.2. Let S = e~ J^o '^i^t) g{e,Xt) discounted asset price. Then the 

following holds: 

(1) Suppose that 9 G y° , F{9) = / and R{9) = A. Then {St)t>o is a true 
martingale under any P^,a; G D. 



■^It follows from Theorem 12. 141 that if some solution exists on a non-empty interval [0, T], so does 
the unique minimal solution. 
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(2) Let X G D° . The process {St)t>o is a true ¥^ -martingale if and only if 
9 G y, F{e) ^ I, R{e) ^ X and (i{t, 0, -1) = and Ht, 0, -1) = are the 
unique minimal solutions of the Riccati equations (|3.2a[) - (|3.2b[) . 

Using p.ip it is clear that 5* is a P^-martingale if and only if Qr.tgix) = g{Xt) 
for all < t < T and with g{x) = e^^'^K Applying Theorem 12.141 to the extended 
process Z the above result follows immediately 

3.3. Moment Explosions. Here we set i = for simplicity. It is well understood 
that the existence of moments E [S^] with y G M"* is intimately connected to the 
shape of the impli ed volatili t y surface derived fro m the prices of options on the 
underlying 5* (cf. iLed (2004 1. iKeller-Ressel |2011 |). Of particular interest is the 



time of moment explosion, i.e. the quantity 

T+(y) = sup {t > : E [S^] < oo} . 
Applying again Theorem 12. 141 we obtain the following 

Proposition 3.3. Let S = cxp {0, X) with and let y G R'' . 

(1) If y + G y° then T^(jj) is the maximal lifetime of the solution {p,q) of 
the extended Riccati system 

(2) If y + G y then T^{y) is the maximal lifetimes of the unique minimal 
solution (p, q) of the extended Riccati system. If y + ^ y then T^{y) = 0. 

Related applications incl ude the approximation of more complicated payoff func- 
tions by 'power payoffs', see lCheridito and Wuealter 2012 and portfolio optimiza- 



tion involving power utility, see Muhle-Karbg |2009{ and the references quoted 
therein. 

3.4. Option Pricing. In general, European option payoffs are non-linear functions 
that do not fall under the setup of the previous subsction. Numerically expensive 
Monte-Carlo simulations may be avoided by the method of Fourier pricing, if the 
characteristic function (or Fourier-Laplace transform) is given in closed form, cf. 
ICarr and Madan 1999| . This is the case for affinc processes, and the key for ap- 



plying Fourier pricing is our Theorem 12.261 on complex exp onential momen ts. We 
provide here an extension of Theorem 10.5 from the book of [Filipovid . [20091 Chap- 



ter 10], which has been written in the context of affinc diffusions, where certain 
simplifications occur (most importantly 3^° = W^). For general affinc processes 
with jumps we have to impose some stronger assumptions and obtain the following 
result. To allow for multi-asset options we consider a generic payoff g : D ^ M. 
depending on all components of the underlying factor process X . In typical appli- 
cations g will be of the more specific form g{x) = ft,(e^*'^^ with h : M^o — » IR which 



can b e accomodated in the theorem below by setting q = 1, see also [Filipovic 
20091 . Theorem 10.6]. 



Theorem 3.4. Let X be an affine process satisfying Assumption \ 2.24\ Assume 
there exists a d x q matrix K such that the payoff function g satisfies 

(3.4) g{x) = f e<-+'^^^->5(A)dA 

for some integrable function g : ^ <C, q < d and with w G 3^°. Suppose that 
(|3.2ap - (j3.2bp has solutions on [0, r] for initial data u — and u — v, which stay in 
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3^° for all t < T. Then we have 



where {(j), ip) are the unique solutions of (|3.2ap - (j3.2bp with complex initial data 
V + iKX. 

4. Proofs for Real Moments of Affine Processes 

4.1. Decomposability and Dependency on the Starting Value. Definition l2.2l 
of an afHne process, imme diately implies a de composability prop erty of the laws 
on the path space (see also lDufhe etall (20031 Thm 2.15]). As in lDufhe et all [2003 . 
Def. 2.14] we write P * P' for the image of P x P' under the measurable mapping 
(w, uj') ^ Lu + Lu' : {n X n,T X ^ {n, T). 

Proposition 4.1. Let X be an affine process with state space D. Its probability 
laws P^ satisfy the following decomposability property: Suppose that x, ^ and x + ^ 
are in D. Then 

(4.1) * P« = P" * P^+^ 

Proof. Write u = {u^, . . . , u^) for an ordered set of points u'^ G 14. Choosing some 
finite sequence < ti < • • • < <Ar in K^O: define 

fix,u)=e^p(j2{Xt,,u'^)^ , (.TG A("\...,«^)ew^), 

i.e. f{x, u) is the joint characteristic function of {Xt-^ , . . . , Xt^) under P^. Applying 
the afhne property (|2.3p recursively we obtain 

(4.2) /(x,u) =cxp(p(u) + (x,<z(u))), 
where p{u) = p\ and (z(u) = qi, with 

(4.3) = </>(tfc +u'=) PAr=0, 

(4.4) =V(tfc-tfc-i,<7fc+w'=), 9^=0. 
From (|4.2p we derive that 

/(a;,u)/(C,u) = /(G,u)/(x + ^,u) 

for all u = (u^ . . . , M^) G . Since the distribution of a stochastic process is de- 
termined by its finite-dimensional marginal distributions, this equality is equivalent 
to gH). □ 

In the following we set 

g(t,y,x)=W [e<^^^'>l = / e<^^«>pt(x, d^, 



for all (i, y) G M^o x and x G D. Note that g(t, y, x) is always strictly positive, 
but might take the value -l-oo. By approximating g(t, y, x) monotonically from be- 
low by bounded functions and using the Chapman-Kolmogorov equation we derive 
that 

(4.5) g{t + s,y,x)= g{t,y,^)ps{x,d^) 

J D 
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holds for all i, s £ K5.0, y G M'' and x ^ D, where +00 is allowed on both sides 
and in the integrand. The following Lemma concerns the role of the starting value 
Xq = X of the affine process with regards to finiteness of exponential moments. 

Lemma 4.2. Let X be an affine process on D and let {T,y) G K^o x K''- Then 
the following holds: 

(a) E° [e^y^^-r^] = 00 implies [e^^'^-^'^^] = 00 for all x G D° . 

(b) W [e^y^^-^)] < 00 for some x e D° implies E"" [e<?^'^^>] < 00 for all x e D. 

(c) W [e'^y^^-r)] < 00 for allx e D implies [e<y^-^'>] < 00 for all t G [0,T],x G 
D. 

Proof. As before we set g{t, y, x) = [e^'^''^'^] which takes values in the extended 
positive half-line (0, c«]. Since D is convex and contains 0, and using the decom- 
posability property of X (cf. Prop. 14. ip we have 

(4.6) g{t,y,Xx)g{t,y,{l~X)x)=g{t,y,0)g{t,y,x), xgAAg[0,1]. 

Clearly g{t, y, 0) = 00 implies g{t, y, Xx) = 00 for all x G I? and A G (0, 1) and we 
have shown (jg). We prove (|b| by contraposition: Assume that g{t,y,x*) ~ 00 for 
some X* G D. We want to show that g{t, y,x) = 00 for all x £ D° . If g{t, y, 0) = 00 
we can use (jaj) and are done. Assume therefore without loss of generality that 
g{t,y,0) < OQ. Equation ()4.6p implies that 

g{t, y, x*/2f = g{t, y, 0)g{t, y, x*) = 00, 

and by iteration it is obtained that g{t, y, 2~''x*) = 00, for all fc G N. Let now x be 
an arbitrary point in D°. Clearly, there exists some large enough fc G N such that 
also ^ = X — 2~^x* is in D° . Using the decomposability property (|4.ip . we obtain 

9{t, y, 0)9it, y, x) = g{t, y, 2~''x*)g{t, y, ^) = 00. 

and (|b| is shown. 

To show (jcj) pick an arbitrary x G D° and e > 0. Since X has cadlag paths 
we can find 5 > such that F^'iWXt - x|| < e) > i for all t < d. With pt{x,d(,) 
denoting the transition kernel of X and B^{x) the open ball of radius e around x 
we can rewrite this as pt{x, Bf{x)) > ^ for all t < 6. We show assertion (jcj) for 
t € [T — i5, r] ; the general case follows then by iteration. By (j4.5p 

9{T,y,x)^ / g{t,y,£)pT-t{x,dCj 

J D 

holds for all {t,y) G [0,r] x R'*. By assumption the left hand side is finite, and we 
want to show that also g{t,y,^) is finite for all ^ G D. Assume for a contradiction 
that 5(t,y,r) = 00 for some t G [0,r] and ^* G D. Then by Lemma W/I{b\ 
git^y^O = 00 for ah ^ G D°. But pT-t{x,D°) > pT-t{x, B^{x)) > \ and we 
conclude that g(T,y,x) ~ 00 which is a contradiction. □ 

4.2. From Moments to Riccati Equations. In this section we prove Theo- 
rem |2TT^1 except for the minimality property of the Riccati solution. 

Lemma 4.3. Let X be an affine process on D, and let T > 0. Suppose that for 
some X e D° and y G R'' it holds that E^ [e<2'"^^>] < 00. Then y e y and the 
following holds: 

(a) There exist functions t iH- p{t,y) G M and t q{t,y) G M'' such that (|2.14p 
holds for allxeD.te [0,r]. 
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(h) W 



,{q{T^t.v),Xt) 



] < 00 for all t G [0, T] and 



(4.7) 



Ft\ = exp (piT -t,y) + {q{T - t, y), X,)) for all x e D. 

(c) The functions p{t,y),q{t,y) satisfy the semi-flow equations 
(4.8a) piT,y)=p{T~t,y)+p{t,q{T^t,y)), p(0,y) = 0, 

(4.8b) q{T,y)^q{t,q{T~t,y)), <z(0, y) - y, 

for all t e [o,r]. 

Proof. From LemmaEait follows that [e^^'-^'^j < oo for all (t.x) e [0,r] x D. 
Fix t e [0,T] and write g{t,y,x) = W [e^^^^^')]. Then by Proposition liH o(t. v. x) 
satisfies the functional equation (|4.ip . Since g{t,y,0) > there exists p{t,y) G K 
such that g{t,y,0) = eP(*'f). Set h{t,y,x) = e'P^^^y'^ g{t,y,Q). Then h{t,y,x) is 
finite for all x £ and satisfies Cauchy's functional equation 

h{t, y, x)h{t, y, £,) = h{t, y,x + £,), x,^,x + ^ £ D. 

We conclude that there exists q{t,y) G M'' such that h{t,y,x) = e^''^*'^^'^^ for all 
X G D, and we have shown (|2.14p . 

To show equation (|4.7p note that by the Markov property of X 



e^y'i'>Pt.,{X,,dO 



holds for all n E N,x E D,0 < s < t. Using dominated convergence we may take 
the limit 7i ^ cxd and obtain equation (|4.7p from (|2.14p . Taking (unconditional) 
expectations in (|2.14p yields 

exp(p(r, y) + {x,q{T, y))) 

= exp {p{T -t,y)+ p{t, q{T - t, x)) + {x, q{t, q{T~t,y)))), 

for alH > and x E D. Since D contains and linearly spans M'', the semi-fiow 
equations (|4.8p follow. □ 



Note that if t H- p{t, y) and t H> q{t, y) are differentiable with derivatives F{y) 
and R{y) at zero, then it follows by differentiating the semi-flow equations (|4.8p 
that (p, q) is a solution of the extended Riccati system (j2.1ip . The main diffi- 
culty thus is showing the differentiability of p and q. This is very similar to 
the regularity problem for affine processes, where the same question is asked re- 
garding the functions 0(i, u) and u) in Defi nition 12.21 Several solutions of 
the regularity problem have been given, see e.g. Keller-Ressel et al.l 2011bl | and 



Keller-Resse 



of Cuchiero 



Schachermaver. a.d Te l^^:^ [2oAj . Here we adapt th7;^;oach 



2011 



to our setting. 



We enlarge the probability space (fi, F, F, P) such that it supports d + 1 inde- 
pendent copies of the affine process X , which we denote by . . . , X''. Without 
loss of generality it can be assumed that X = X^ . In what follows we will use the 
convention that upper indices correspond to the different instances of the process 
X, while lower indices correspond to the coordinate projections of a single pro- 
cess. For a vector G denote by the restriction of the probability law P 
to {X\Q) =x'], i.e. the process X^ starts at the point with P^ -probability 1. 
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Similarly for an ordered set x = . . . x'') of points in D, we denote by P'' the re- 
striction of P to {(X°(0) = xO) A • • • A (X'^iO) = x'^)}, i.e. the processes . .X'* 
start at the points x'^,. . .x'^ respectively with P'^-probability 1. 

Lemma 4.4. Let X^, . . .X"^ be d + 1 independent copies of the affine process X . 
Furthermore let x = (x*^, . . . x'') he d -\- 1 affinely independent points in D . Define 
the matrix-valued random function 



^(t;x,w) = 



1 XUt,oj) 



1 Xi{t,u:) 




Then there exists 5 > Q such that 

P^(dctS(t;x) 7^ for all < t < 5) > ^, 

i.e. t S(t;x) stays regular on [0,(5] with ¥^ -probability at least i. 

Proof. Since {x'^,...x'^) are affinely independent, and X*(0) = x* for all i £ 
{Q,...d} with P'^-probability one, the matrix S(0;x) is regular P'^-almost surely. 
Define at = infse[o.f] |dctS(s;x)|. Since the processes X' arc right-continuous, also 
det S(s; x) is, and hence even a^. By dominated convergence also bt = V^{at > 0) is 
right-continuous and has the starting value b^ = 1. We conclude that there exists 
some S > such that bs > i, which completes the proof. □ 



The following proposition settles Theorem I2.14[|5| apart from the minimality 
property of (p, q) as solutions of t he extended Ric cati system. The key ideas in the 
subsequent proof are coming from lCuchiero 2011 . Proofs of Lem. 1.5.3, Thm. 1.5.4]. 



Proposition 4.5. Let X be an affine process on D and let T > 0. Let y S M'' and 

suppose that [e^^'''^'^^] < oo for some x G D° . Then y £ y and there exist a 
solution {p,q) up to time T of the extended Riccati system (|2.11[) . such that (|2.14p 
holds for all x £ D, t e [0,T]. 

Proof. Recall that we are working on an extended probability space that supports 
d-\-l independent copies {X^, . . . , X"^) of X. Let x = (x°, . . . , x'') be c? -I- 1 affinely 
independent points in D. By Theorem 12.41 each X' is a P'^-scmi- martingale with 
canonical semimartingale representation 



(4.9) 



XI ^ x ' 



+ / b{Xl_)ds + Nl + 



i^~h{0)J'{u;;ds,dO, 



where is a local martingale and J''{uj;dt,d^) is the Poisson random measure 
associated to the jumps of X' with predictable compensator K{Xl_, d^)dx. 



By Lemma 14.31 we know that 



(4.10) 



cxp{piT~t,y) + {q{T~t,y),Xi)) 



for each i e {0, . . . ,d} ,t e [0,T]. Let us denote = E^ 



Clearly, 



each t ^ Ml is a P^-martingalc for t < T and for each i G {0, . . . , d}. Taking 
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logarithms and arranging the equations in matrix-form we get 

/ P{T-t,y) \ 
qi{T - t,y) 

and recognize on the right hand side the matrix S(i;x, w) from Lemma 14.41 The 
latter allows to conclude that there exists a set A (Z VI with W'^{A) > i and some 
S > such that S(t;x, w) is invertible for all t G [0,(5] and uj G A. Hence for 
T' = TA5we obtain 
(4.12) 

\ 1 Xf{t,u) ■■■ Xi{t,u) 

for all t G [0,r']. All processes occurring on the left-hand side of equation (j4.12p 
are semimartingales, hence also the right-hand side consists row-by-row of semi- 
martingales for all t S [0, T']. Since they are deterministic, the functions t p{t, y) 
and 1 1— > q{t, y) are of finite variation on [0, T']. This implies in particular that they 
are almost everywhere differentiable and can be written as 

(4.13a) p{T' -t,y)- p{T' , y) = - f dp{T' - s,y), 

Jo 

(4.13b) q{r -t,y)- q{T', y) = - ( dq{T' - s,y), 

Jq 

Applying Ito's formula to the martingales M^'^ we obtain 

Mi'y = M^'y + f Ml^l {-dp{r -s,y) + {-dqir - s,y),Xl_)) + 
Jo 

+ Ml-y ( {qir - s, y), b{Xl_) + \ {q{T' - s, y), a{Xl_)q{r - s, y)) + 
+ ^ (e('(^'-^^^)-«> - 1 - {q{T' - ,5, y), /i(C))) K{Xl_,dO)ds+ 

+ f Ml'y{q{T'-s,y),dNl) + 
Jo 

^lo Id^^^'- (^^'^^'"''"''^^ - 1 - - {J{u^,ds,dO-K{Xl_,dOds) , 

for alH G {0, . . . , d}. On the right hand side, the last two terms are local martingales 
and the other terms are of finite variation. Hence the finite variation terms have to 
sum up to 0. Rewriting in terms of the functions F(y) and R{y) this means that 

-dp{T'-t,y) + {-dqiT' -t,y),Xl_) = FiqiT'-t,y))dt + {Xl_,RiqiT' -t,y)))dt 



(4.11) 



/ logAf,"(a 



V logA/,'*(c^) 





' log MO (c.) 
logA/t''(a;) 



/ v{T'-i.y) \ 

qi{T'-t,y) 
V qd{T'-t,y) J 
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holds for almost all t £ [0,T'] P'^-a.s. Inserting into (|4.13|) and using the regularity 
of the matrix S(t,x, w) on [0, T'] this yields 



(4.14) 
(4.15) 



piT' -t,y)-p{T',y) = - f F{q{T' 

Jo 



,y))ds, 



q{T' ~t,y)~q{T',y) 



R{q{T' - s,y))ds. 



Applying the fundamental theorem of calculus, we have shown that (p, q) is a so- 
lution to the extended Riccati system (|2.1ip up to T' = T A S, where S was given 
by Lemma l4.4l To show the general case wc conclude with an induction argument. 
Suppose that {p, q) are solutions of the extended Riccati system up to Tk ~ TA{k6). 
We show that they can be extended to solutions up to T^^i = T A {{k + 1)S). Set 
Afe = Tfc+i - Tk; clearly Afe < S. By Lemma [e<2'-^^>] < oo implies that 



< cxD and by Lemma 14.31 we have that E 



< oo. Set 



y' — liy^ Tk). Then - proceeding exactly as in the proof above - we obtain 



(4.16a) 
(4.16b) 

for all t e [0, Afc] 

(4.17a) 
(4.17b) 



-p{t,y') = F{q{t,y')), 
§-^q{t,y') ^ R{q{t,y')), 



p{0,y') = 
9(0,2/0 = 2/' 



Using the flow property this is equivalent to 
d 



dt 
d_ 

dt 



pit,y)) = F{q{t,y)), p{0,y) = 
q{t,y))^R{q{t,y)), q(0, i/) = 2/ 



for all t e [Tk,Tk+i]. By the induction hypothesis (j4.17p already holds for all 
t € [0, Tk] and we have shown that (p, q) is a solution of the extended Riccati 
system up to Tk+i = T A S{k + 1). As this holds true for all fc G N the proof is 
complete. □ 



4.3. Prom Riccati Equations to Moments. Using the result from above, the 
step from the extended Riccati system to the existence of moments is simple: 

Proposition 4.6. Let X be an ajfine process taking values in D. Let y £ y and 

suppose that the extended Riccati system (j2.1ip has a solution {p, q) that starts at 
y and exists up to T > 0. Then W [e<^^-^'^>] < oo and ([2ll| holds for all x E D, 
t £ [0,1"], where {p,q) is also a solution up to T to (|2.1ip . 

Proof. Using the solution {p, q) of the extended Riccati System (|2.1ip define for 
tG [0,T], 



(4.18) 



= exp {p{T ~t,y) + {q{T - t, y),Xt)) . 
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Applying Ito's formula to Mj^ and using the semimartingale representation ()4.9|) 
we see that 



Ajy + I Ml_ {~p{T -s,y) + {-q{T - s, y),Xs^)) ds+ 



M. 



, y), biXs-)) + I {q{T - s, y), a(X,_)g(T - s, y)) + 



+ / (e<«(^-«.i').?)_l-(g(T-s,y),MO))i^(X,_,dC))ds+ 



Ml_{q{T~s,y),dN,)- 



1 - {q{T - s, y), h{0) {J{^. ds, dO ~ KiX,^,dOds) . 



The ds-tcrms can be simplified to 

-p(T ~s,y) + {-q{T - s, y), + F(^(T - s, y)) + {R{q{T - s, y)), - 

and we conclude that (Mf )(g[o^T] is a local P^-martingale for all x E D. It is also 
strictly positive, and hence it is a P^-supermartingale. Therefore 



< M-j < oo 



for all X G D. The second part of the assertion, and in particular the validity of 
equation (|2.14p follows now by applying Proposition 14.51 □ 



4.4. Proof of Theorem 12.141 Looking at Proposition 14.61 and Proposition 
we see that Theorem 12.141 is almost proved. Only one issue in both parts of the 
theorem is not answered yet, namely the minimality of (p, q) in p.l4p as minimal 
(hence unique, see Remark l2.13p solution of the extended Riccati system. We start 
with the following lemma: 

Lemma 4.7. Let {p,q) and {p,q) be given as in Proposition \4.6\ Then for all 
t e [0, T] andx e D 

P{t, y) + y), x) < pit, y) + {q{t, y),x) 

Proof. Set Mf = exp {p{T - t,y) + {q{T - t,y),x)) and define Mf as in jpi^ . 
Then, for each x ^ D the process is a P^-martingale (see (|4.10p f): is a 
P^-supcrmartingale, and they satisfy A/^ = M^. Hence 



Mf 



w[My\Tt] 



for all t G [0,T]. Taking logarithms the claimed inequality follows. 



□ 



Proof of Theorem \2.14\ Proof of (jlj) : In view of Remark 12.131 we only need to 
show that the solution {p, q) of the Riccati system established in Proposition 14.51 is 
minimal. Let (j},q) be another solution on [0,T'] of the extended Riccati system, 
T' < T. Then by Proposition there exists {p*,q*) such that ^^TT^ holds for all 
y G D and t £ [0, T'], as is the case for (p, q). By taking logarithms of the respective 
right sides of (|2.14p and by applying Lemma l4Jl we see that 

p+ {Q,y) =p* + {q*,y) <p + {q,x), 
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on [0,r'] and for all y E D. Hence by Definition 12. Ill (p. q) is the minimal solution 
of the extended Riccati system, and we are done with part (jaj). 

The proof of (|b| follows immediately from Lemma I4.7[ Definition 12.111 and Re- 
mark [231 

□ 



5. Proofs for Complex Moments of Affine Processes 



In this section wc show Theorem 12.261 on the existence of complex moments 
of affine processes, whose state space satisfies Assumption 12.241 The key to the 
proof is to relate the life time of solutions (0, ip) of the complex Riccati system 
(|2.15ap - (|2.15b[) and the solutions {p, q) of the extended Riccati system (|2.11ap - 
(|2.11b|) . Unlike in preceding parts of the paper, we only solve for initial values in 
the interiors y G y° resp. u € S{y°). Also, in this section we need more precise 
knowledge about the restrictions on the parameters, which appear in the Riccati 
equations. 

With iSj" we denote the d x d positive semidefinite matrices. Let T+{y) resp. 
T-|_(u) denote the maximal lifetime of i i~> y), q(t, y)) resp. t {(t>{t, u), -tplt, u)). 



Proposition 5.1. Suppose that Assumption \2.24] holds true, and let u £ S{y°) 
and y — Re(u). Then T^{u) > T^{y). 

We split the proof into the two cases covered by Assumption 12 . 24l ~ a state space 



D of the form 



and a state space of the form . Note that 



St are convex cones. To apply certain results of IVolkmannI |1973l | on multivariate 



and 



ODE comparison, wc introduce the following property: 

Definition 5.2. Let K C M'^ be a proper closed, convex cone, and denote by ^ the 
induced partial order. Let U C M''. A function f : U — > R*^ is called quasimonotone 
increasing (with respect to K), if for all y,z £ U for which y ^ z and (y, x) = (z, x) 
for some x € K it holds that {f{y), x) < {f{z), x). 



5.1. State Space D 



^^0 



In this section we consider the 'canonical 
We use the index sets 



from Duffie et al. 1200. 



state space' D 

I = {1,2,..., m} and J = {m + 1, . . . , d} corresponding to the positive and to 
the real valued components of D respectively. Accordingly, Ri denotes the func- 
tion . . . , and similarly Rj is constituted by the last n coordinates of 
R. 

First, we recall the definition of the admissible par ameter set for (conservative) 
affine processes on 



from Duffie et al. 2003 



Definition 5.3. A set of ]R''-vectors b, P^, . . . , (3'^, positive semidefinite d x d ma- 



trices a, a^, . . 



D = 



iff 



Levy measures m, fi^ , . . . , fx"^ on M'', is called admissible for 

aki = for all fc G / or / G /, 
a-' ^ for all j G J, 
all = if fc G / \ {i} or / G / \ {i}. 
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beD, 

J > for alH G /, fc G / \ {i}, 

PI = for all j eJ.kel. 



PI 



I |^/|m(d^) < oo, 

fjL^ = for all j e J, 



'kl<l 

Remark 5.4. The matrices a,a^ , . . . , a"^ are frequently referred to as diffusion matri- 
ces, the vectors &, . . . , as as drift vectors and the Levy measures m, fi^ , . . . , fi''' 
as jump measures. 

Let R{y) = {Ri{y), . . . , Rd{y)) be defined as in Proposition [2]8] The admissibil- 
ity conditions imply that each Ri(y), . . . , Rd{y) is a convex lower semicontinuous 
function of Levy Khintchine type. Denoting ^^[d^f] :~ in[d(f) we therefore have 

(5.1) 3^= jtyeM'*: V / e<^^«VXrfO < oo| 
which is the intersection of the effective domains oi F, Ri, . . . , Rd- 

We start with the following crucial lemma: 

Lemma 5.5. There exists a function g which is finite, non-negative and convex on 
y such that for all u G >5'(3^°) we have 

(5.2) Re {{uj, Ri{u))) < g(Re u){l + \\uj\f){l + \\uif). 
Proof. It clearly sufficient to show 

Re {u^R^{u)) < g,(Re(u))(l + |w,/P)(l + lu/P) 

individually for each j G / and with some non-negative convex gi{.) that is finite on 
y. In addition we may split Ri{u) into the drift part, the diffusion part, a small- 
jump part and a large-jump part and show the inequality for each part separately. 
The drift and the large jump-part are the easiest to deal with. Using the Cauchy- 
Schwarz inequality we infer the existence of a positive constant C such that 

(5.3) Re{M{P\u))) < \u,\m\\^^i\ + W'j\\ui\) < 

<C(l + \uj\^)(l + \uj\') 

and 

(5.4) Re lu, [ e<^'"^^l\do] < \u^\ [ e<«'««"V*(rfe) < 

<5,(Reu)(l + |u/n 

for the large-jump part. Here gi{z) := e^^'^V*('^C) clearly is a non-negative 

convex function which is finite on y. To estimate the diffusion part we have to take 
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into account the admissibility conditions, which tell us that alj is zero if j G /\ {i}. 
Thus we obtain 

(5.5) Rc (ui (u, a'w)) = a^Juip Re + 2 Re (juipaijuj) + Re (uiuJajjUj) < 

<C(l + (Rc ui)+){l + \uj\^){l + \ui\^), 

as desired. The harde st term to estimate is the small-jump part. We follow the 
proof of Lemma 6.2 in iDufBe et al.l j2003j . As a shorthand notation we introduce 
u/_ — ui\{i} and uj^ = wjuji}- First we do a convenient Taylor expansion of the 
integrand h{£,) ^ e<«'"> - 1 - {ij+,uj+) with |^| < 1: 

(5.6) h{0 = e<«'"> - e<«'^+'"^+> + e«*"' (^e<«-^'"-^> - 1 - (C/, uj)) + 

+ (O, «j) (e«-"' - 1) + - 1 - ^,u, = 

Next we calculate 

Rc {mh{0) = K{u,0 + [ (1 - i) Re (u,e»e*"'«0 



Since |^| < 1, we get 

(5.7) \K{u,0\ < e^""" d^^.!' + \ui\\uj\' + \ui\\uj\) + < 

< (1 + e(«-")+)(l + \uj\'){l + \uj\') (16-1 + n 
for the first term. For the second term we use Lemma 15.61 below and estimate 



(5.8) |u,f J (l- t)Re (u,;e;e*"'«') dt < \u,\^^, (e«'(«°"')+ - l) . 

Adding up (|5.7p and (|5.7I) . and integrating against the Levy measure /z' we obtain 

(5.9) Re (^^J^^^^ h{^)fi'{dO^ <UReu)il + \uj\^){l + \ui\'') 
with 

j\i\<i J\^\<i 

which is non- negative, convex and finite for all y 6 R'^. Adding up (|5.3p . (|5.4p . 
((53)) and ((O)) yields the desired estimate ((5?2|) □ 



Lemma 5.6. For any z G 



(5.10) (1 -<)Rc(ze*^)rfi < (e(^°")+ - l) 
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Proof. For Re z < the inequality was shown in iDuffie et alJ 2003l | . Denote the 
left hand side by L{z). Writing z = p + iq and evaluating the integral we have that 



L{z) 



(1 - f)e^'* {pcos{qt) -qsm{qt)) 



1 



{p (eP cos{q) -1 - p) + q (e^ sin(g) ~ q)} ■ 



The expression is symmetric in q such that we may assume that q > 0. If in addition 
p < then using cos{q) > 1 — q'^/2 and sin((7) < 1 we may estimate 

Liz) < {p{eP ~l^p)-q' + q^eP{l - p/2)) . 

Since e^(l — p/2) < 1 and (e^ — 1 — p) > the right hand side is smaller than 
showing the Lemma for p < 0. li p > we may use that cos{q) < 1, sin(q) < q and 
eP — 1 < peP to estimate 

Liz) < {p'ieP - 1) + q\eP - 1)} = ie^ - 1), 



□ 



p^ + q 

thus completing the proof. 

Recall Definition 15.21 of quasimonotonicity with respect to a convex cone K. 
Here, K = R™; in this particular setting, quasimonotonicity of a function f : U C 
K — > R™ can be expressed in coordinates and is equivalent to : 

y ^ z, and yi = z^ for some z G {1, . . . , m} ^ /^(y) = 

Lemma 5.7. Let yj € M". For each t > 0, yi t-^ Riiyi,ipjit,yj)) is quasimono- 
tone increasing (with respect to the natural cone W^) on y. 

Proof. See, for instance, Keller-Ressel )2009| or Maverhofer et al. 2011 

We further need the following special property of y° : 

Lemma 5.8. If y E y° , z G M"^ and zj ^ yj, zj ~ yj, then we also have z ^ y° . 

Proof. We choose e>0 such that B^^^iv) ^ [w \ \\y - w\\ < e} C 3^. By ([51 
we have for z = 0, 1, . . . , d 

< oo, \\w\\ < £■ 



□ 



(5.11) 



l«l>i 



Note the Levy measures /i* are clearly positive and supported on D. Now for all 
^ G D we have 

(z + w, = zj^i + zJCj + (w, = z]^i + 2/Je/ + (w, < + 

because ^/ G R"' and Zi < yi for all i G /, by assumption. Hence, by the mono- 
tonicity of the exponential wc sec that (|5.1ip holds with y replaced by z. Hence, 
once again by (|5.ip we have i3<e(z) C y, i.e. z G 3^°. □ 

We are now prepared to prove Proposition 15.11 under Assumption 12.241 (a): 

Proof of Proposition \5. 1] under Assumption \2.24\ (a). By a straightforward check, 
for every u € y, 

Re(i?i(u)) < i?i(Re(u)), 

and by Lemma ISTfl we can apply the ODE comparison result of lVolkmannI 1973l | to 
the first m coordinates of ijj, which let us conclude that Re(V'/(t, u)) ^ ipiit, Re(M)) 
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for t < T+{u) A r+(Re(u)). In view of Lemma 15.81 we therefore have T+(u) > 
r+(Re(M)), unless explodes before Re(M))|| does. We show in the 

following that this cannot happen: By Lemma 15.51 we have 

< g(Re u)){l + \\^Pj{t, u)f){l + ||^7(t, u)f) 

with a function g which is finite on all of y. Since ipj{t, u) = ipj{t, uj) exists glob- 
ally as solution of a linear ordinary differential equation, we obtain by Gronwall's 
inequality applied to (1 + ||'07(^, u) ||^) that 

\mt.u)\\ < \\uir + {i + \\ui\\') 

where h{t) := g{Re ^|;{t,u)){l + \\tl;j{t,u)\\'^). Hence we have shown T+(u) > 
r+(Re(M)). □ 



Hence we have shown T^{u) 



5.2. Matrix State Spaces. Let 5*^ be the space of symmetric real dx d matrices, 
endowed with the inner product (x, y) = tr(xy), where tr denotes the trace operator. 
We further denote by C™^" the space of complex m x n matrices. We make the 
latter into a normed space by introducing a norm as ||a||^ := tr (aa^). Here ^ 
denotes matrix transposition, and a is the element-wise conjugate of t he matrix a . 
W e start with the following observation, which is a generalization of Maverhofer , 



20111 . Lemma B.l]: 



Lemma 5.9. 

all a e C'"^" 

(5.12) 



There exists a locally Lipschitz function h : Sd 
and for any b £ S{Sn) we have 

Rctr(-6aa^) < h{Rcb) ■ ||af . 



such that for 



Proof. Recall that the projection tt : S'rf — > S'J is a well defined, convex (hence 
locally Lipschitz continuous) map, which satisfies 7r(z) >: z for all z S §(j. 

Let us write a = a\~\~ ia2 and b = bi + ib2 with ai, 02 G jj'^x" and 61, 62 G Sn- 
Then we have 

Retr(— 6a^a) = Rctr(— (61 + ib2)iaj — iaj){ai + 102)) 

= ti{-bi{aj ai)) + tT{-bi{aJ 02)) + 

< tr(7r(-6i)(a7ai)) + tr(7r(-6i)(aja2)) 

<lkMi)ll(l|aif + li«2|P) 

^M-bi)\\ ||af. 

The last but one inequality holds in view of the Cauchy-Schwarz inequality, 
now see that inequality (|5.12p holds by setting h{x) := ||7r(— a::)||. 



We 

□ 



Next we present the admis sibility conditions for matrix- valued affine processes 



201lj . Note that in the case d ~ 1 



that have been established in iCuchiero et al. 

it holds that S'^ = Mjjo, i.e. the one-dimensional case is already covered by the 
previous section. Therefore we may assume that d > 2, which le ads to several sim- 
plifications of the parameter conditions. It has been shown in iMaverhofen 2011 
that affine processes on St (d > 2) do not exhibit jumps of infinite total vari- 
ation. Compared with ICuchiero et al.l [201l| this makes the use of a truncation 
function in the definition of R obsolete and also simplifies the very complicated 
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( because hard-to-check ) necessary tradeoff between linear jump coefficient and drift 
( (Cuchiero et all [20111 (2.11)]). In the following, ^ denotes the partial order on Sd 
induced by the cone S'j'. 

Definition 5.10. An admissible parameter set {a, b, B , in{d(,) , fi{d£^)) consists of 

• a linear diffusion coefficient a G S'^ , 

• a constant drift b £ S'^ satisfying 

by {d- l)a, 

• a constant jump term: a Borel measure m on S'^ \ {0} satisfying 



A l)m{d^) < oo 

s,+ \{0} 

a linear jump coefficient jj, which is an S'J"-valued, sigma-finite measure on 
\ {0} satisfying 



A l)n{d^) < oo, 

• and, finally, a linear drift B, which is a linear map from Sd to Sd and 
"inward pointing" at the boundary of S'^ . That is, 

tT{xB{u)) > for aU u,x e Sj with tr(w.T) ^ 0. 

Remark 5.11. Using the notation a(x) from p.7p we have a{x){u) ~ 2tT(xuau), 
and the following are equivalent: 

(1) Condition [121 dnl) 

(2) a = or a is invertible. 

: all X G 

^d \{0}- 

The only non-trivial direction to prove is ([T]) =^ ([2]). Assume, for a contradiction, 
that a 7^ 0, but a is degenerate. Then there exists u E \ {0} such that ua = 
au = 0. But then a{x){u) = ti{xuau) = 0, for any x. 

Note that Cuchiero et al.l 2011 1 uses the Laplace transform to define the affine 



(3) Either a{x) vanishes for all x € S'^, or it is non-degenerate for any x £ 



property, which introduces seve ral changes of signs c ompared with our definition. 
To comply with the notation of Cuchiero et al.l |201lj we introduce 



F{y) = -Fi^y), R{y) - -R{-y), 
which can now be written as 

Fiy) = triby) - / _ l) m(dC), 

Jsi\{0} ^ ' 

R{y) = -2yay + B^ (y) - [ (e- - l) ^i{dO 

Js+\{o} ^ ' 

Writing furthermore 

V) = -P{t, -y), Q{t, y) = -q{t, -y), 
and similarly for (/) and t/;, then, bv lCuchiero et al.l j201l| . 



g-tr(yXt 



g-p(t.J/)-tf(9(t>J/)2:) 
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for all i > 0, y, X e and by iMaverhofed 201l| the exponents (p, q) : R+ x 5 -> 
R+ X Sj" solve the system of generalized Riccati equations 

(5.13a) ^^p{t,y) = F{qit,y)), 

(5.13b) ^q{t,y))^R{q{t,y)), 

given initial data p(0, y) = 0, g(0, y) = 

Since /i is an S'j'-valued measure, tr(/i) is a well defined nonnegative measure, 
naturally given by 

tr(A.)(A)=tr(/i(A)). 
Accordingly, the domain y :— —y is given by 

(5.14) y^\yeSd\f e-'<y^\m{dO + tr{li){dO) <oo\ . 

The inclusion D holds in view of the positiv e definiteness of the measure ^, while the 
inclusion C follows from [Maverhofer , 2011 . Lemma 3.3]. Similarly to the preceding 
section, we start with the following crucial estimate. / denotes the dxd unit matrix. 

Lemma 5.12. Suppose that the diffusion coefficient satisfies a ~ I or a = 0. 
Then there exists a locally Lipschitz continuous function g on y° such that for all 
u € S{y°) we have 

(5.15) Re(tr(u^(it))) < g(Re(u))(l + ||uf ). 



Proof. As in the proof of Lcmma [5.5l wc start with drift and big-jump parts. Clearly 
we have 

(5.16) Rctr(t2B^(w)) < Gi(l + ||ii||)2 

for some positive constant Gi. What concerns the big-jump parts, we have 



Retr u 



^ - l) ^,{d0] j < lkl!tr(/i)({C : H^ll > 1}) 



(5.17) +\\u\\ (e-'^(^^<AtT{f^){dO 

Ju\\>i ^ ' 

(5.18) <,92(Re(u)(l + l|u||'), 

for some locally Lipschitz continuous function g^. The integral (|5.17p is finite, 



because by Re(u) G 3^ by assumption. Here we have also used (Maverhoferl 12011 
Lemma 3.3]. Note that we can set 

g^ tr(M)({C : ||e|| > 1}) + / (e" X^A^i)- 

For a = we set 33 = 0. If a = /, we involve Lemma 15.91 and obtain 
(5.19) Re(Mw2) <g3(Re(^))||y||2, 

where gz{-) = h{-) = 7r(--). 
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It remains to estimate the small-jump part. Using again Maverhofen . 12011 
Lemma 3.3] we have 



Retr u 



'0<||«||<1 



(^e-'>^«) - l) ^(de) j = Rctr [s /" / trK)e-" ds^(de) J 




'Il5l|tr(/x)(d^) 

<ell / ||ei|tr(/.)(dO 

Jo<||5||<i 

< 54(Re(w))(l + ) 



with 



94{y):=eM ||^||tr(/i)(dO. 

Jo<||c||<i 

Summarizing the last estimate together with (|5.16p . (|5.18p and (|5.19p and setting 

g{y) := d + 52(2/) + gsiv) + Qiiy) 
proves the assertion. □ 

We provide two further lemmas: 

Lemma 5.13. If y E y° , and z £ such that z ^y, then we also have z e 3^°. 

Proof. Using (j5.14p we infer the existence of some £ > such that for all w € 
B<e{0) = {w e Sd \ \\uj\\ < e} we have 

e-i'^iiy+^K)(m{dO +tT{^l){dO) < 00. 

I«ll>i 

The assumption of the lemma implies that (z,^) > {y,£,) for all ^ £ S'^ . Further- 
more, m and tr(/x) are supported on S'^ . Therefore we have 

e-"-«"+"')«)(m(de)+tr(Ai)(dO) < / "■«y+"')«)(m(d^)-Ftr(/i)(dO) < 00 

IICII>i "'ll?ll>i 

for all w S i?<e(0) which in view of (|5.14p proves that z G y°. □ 
Lemma 5.14. R is quasimonotone increasing (with respect to ) on y° . 

Proof. The proof is analogous to the one of jCuchiero et al. . 20 111 Lemma 5.1], 
which states quasimonotonicity of R on S'^ . □ 

We are now prepared to prove Proposition 15 . 1 1 for D = S'^ , d > 2: 

Proof. According to Cuchiero et al.l . l201ll Theorem 4.14] for any affine process X 
(with diffusion coefficient a) there exists a linear automorphism g of S'^ such that 
the affine process Y = g{X) has diffusion coefficient a = diag(Ir, 0), where is the 
r X r unit matrix, and r = rank(a). According to our assumption r = or r = d 
(see Remark 15. lip , and linear transformations do not affect the blow up relation 
(between the real and complex-valued solutions) we are about to prove here. Hence 
we may without loss of generality assume that a = or a = /. 

For any u S S{y°) we write y = Re(u). The q uasim onotonicity of R (Lemma 
I5.14P let us apply the multivariate comparison result by IVolkmannI |l973j and we 
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conclude that for t < T+{u) AT+{y) we have Keip{t, u) >: q{t, y). In view of Lemma 
I5.13l we only need to show that t h- > ||'!/'(i7 does not explode before 1 1— ?■ \\q{t, y)\\. 
By Lemma |5.12[ there exists a continuous function g such that 

Retr(?2^(u)) < g{Re{u)){l + \\uf), u G S{y°). 

Hence, we have for all t < T+{u) A T+{y) 

^^[\\i,{t,u)f) = 2Retr(V^(t,u)i?(7A(<,u))) < ff(Re(V^(<, u)))(l + \m,u)\\^) 
and by Gronwall's inequality, we obtain 

U{t,u)\\<{l + \\uf) f g{s)e^o3{0di^g_ 
Jo 

Hence we have shown that T-^-{u) > T^{y). □ 

5.3. Proof of Theorem I2.26[ The first part of Theorem 12.261 has been subject 
of Proposition 15.11 For the proof of the second part, the validity of the complex 
transform formula (|2.16p . we utilize the concept of analytic continuation. 

Proof of Theorem[KM Fix x e D and t < T. 

Wc proceed differently, depending on the choice of state space. We start with 
D = Sj: Consider the set U := {y e y° \ T+{y) > T}. By the standard existence 
and uniqueness theorem for ODEs, we have that U is an open non-empty set. Also , 
by Lemma r5.14l and Lemma [5.13l and [Keller-Ressel. Maverhofer. and Smirnov 1201(1 



Theorem 1] we have that [/ is a convex, hence connected set. 

Now set U' := S{U) C C^. By Proposition O we have T+(u') > T for all 
u' e U' . Fu rthermore, since u i— >■ R{u) and u F(u) are complex analytic on y°, 
we have by [Dieudonnl . 19691 Theorem 10.8.2] that the function 



M{u) 



e 



— (p{t.u)~tr{ip{t.u)x 



is complex analytic on U'. Furthermore, by Theorem 12 . 141 and by Remark l2.23l we 
have that 

We conclude that the function ^(u) : U' ^ C : u ^ ]Ea:jg-tr(uJft)j ^^^^ analytic 
fun ction, which coinc ides with u M- M(u) on a set of uniquenes^, namely U. Hence 
by iDieudonnel . Il969i (9.4.4)] E=^[e-*' ("^')] = M{u) on all of U' , and we are done 



with the matrix case. 

Let now D = M™ x M". Let q* e y° such that T+{q*) > T. We introduce the 
two sets 

U:^{qey°\ qj = q}, T+{q) > T}, U' = S{U). 
First, we show that the non-empty open set U is convex. For y & U we have 

y) = Rit, X{t, y)), where Rit, •) := (•, ^Pj{t, y))) 
subject to x(0j2/) = Ui- Similarl y to the matrix case the p roperties of Lemma 15.71 



3_pr 

and Lemma [5781 imply in view of Keller- Ressel et al.l . l2010l Theorem 1] that ?7 is a 
convex set. 



this set must be open and connected 



EXPONENTIAL MOMENTS OF AFFINE PROCESSES 



29 



Furthermore for any u G U' , x{t,u) := ipi(t,u) satisfies quite similarly the non- 
autonomous ODE 

^X(i, u) = R{t, x{t, u)), where Rit, •) := Riit, (■, ^Pj{t, u))) 

subject to x(0, m) ~ uj. Since u i— >■ R{u) and u i— >■ F[u) are complex analytic on 
y°, we have that the function 

M{u) g'^'(^") + ('/'(t."),a:> 

is complex analytic on U' . Furthermore, by Theorem 12. 141 we have that 

M{y) = eP(*^^')+<9(*'2')'^) = E^[e<?^-^*>], y € U. 

We conclude now exactly as in the matrix case that the function : U' 

C : u I— > £^[6^"'^^'^] is an analytic function, which coincid es (see Remarkl2.23p with 
u H> Ad{u) on the open and connected set U. Hence by Dieudonne . 19691 (9.4.4')] 
]Ex[g(«,Xt)j ^ ^.^(^j-, jj/^ ^Yie proof is finished. □ 
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